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Abstract. The Heisenberg evolution of a given unitary operator corresponds 
classically to a fixed canonical transformation that is viewed through a moving 
coordinate system. The operators that form the bases of the Weyl representation and 
its Fourier transform, the chord representation, are, respectively, unitary reflection 
and translation operators. Thus, the general semiclassical study of unitary operators 
allows us to propagate arbitrary operators, including density operators, i.e. the 
Wigner function. The various propagation kernels are different representations of 
the superoperators which act on the space of operators of a closed quantum system. 
We here present the mixed semiclassical propagator, that takes translation chords 
to reflection centres, or vice versa. In contrast to the centre-centre propagator that 
directly evolves Wigner functions, it is guaranteed to be caustic free, having a simple 
WKB-like universal form for a finite time, whatever the number of degrees of freedom. 
Special attention is given to the near-classical region of small chords, since this 
dominates the averages of observables evaluated through the Wigner function. 
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1. Introduction 

The semiclassical theory for the evolution of the states of a closed quantum system has 
a long history The way in which operators propagate semiclassically driven by the 
Heisenberg equation is a more delicate matter. On the one hand, observables that give 
rise to smooth, approximately classical functions in the Weyl representation (admissible 
operators pQ) also propagate in an almost classical manner. In contrast, unitary 
operators exhibit semiclassically narrow oscillations in the Weyl representation and so 
do density operators, represented by the Wigner function. The propagation of these 
oscillatory regions is not trivially related to the corresponding classical evolution, but it 
has been shown that the semiclassical Wigner function is transported by pairs of phase 
space trajectories j^j- Alternatively, these may be combined into a single trajectory in 
double phase space, secondary phase space [3 , or, in full rigor, the symplectic groupoid 

A general qualitative picture in terms of a pair of orbits can always be invoked. 
If the chord that connects the initial points of the relevant pair of orbits is small, 
we may substitute these by a single central orbit, so that the propagator for Wigner 
functions assumes the limiting 5-function form along the classical trajectory, deduced 
semiclassically by Marinov |5.. But, in the case of long chords, both the orbits at the 
chord tips must be used. The dificulty is that the initial chord depends on the initial 
state, so that it is not so easy to derive a general semiclassical propagator. Furthermore, 
the limit of small chords is a caustic of the semiclassical theory demanding a higher 
uniform approximation [S| than that developed in E] . 

The Weyl, or centre representation decomposes arbitrary operators into a 
superposition of reflection operators [Zj. The Fourier transform of the Weyl symbol, 
i. e. the chord symbol, is the expansion coefficient of the same operator in the basis of 
translation operators. Both these bases belong to the general class of unitary operators, 
so that the semiclassical theory for evolution of unitary operators allows us to propagate 
arbitrary Weyl symbols, or Wigner functions. The purpose of this paper is to provide 
this general semiclassical framework. If one considers that the space of linear operators 
that act on Hilbert space also forms a linear space, then the kernels of the integral 
representation for the propagation of Weyl symbols and chord symbols correspond to 
representations of super- operators jH]. Here, we are only concerned with a subclass of 
super-operators, or non-selective operations P , appropriate to closed quantum systems, 
but a generalization of the present semiclassical theory to open Markovian systems will 
soon follow. 

The natural setting for the present semiclassical theory is double phase space, 
whose elements are all the ordered pairs of phase space points and hence encompasses 
all possible classical transitions. A uniform translation of the ordinary phase space 
is represented by a plane in double phase space, which is transverse to the double 
phase space plane that defines a canonical reflection through a phase space point. The 
points on these alternative sets of planes can be used as conjugate coordinates for 
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double phase space, which correspond, respectively, to the chord and the centre (Weyl) 
representations. It is the classical evolution of these planes in double phase space that 
supplies the semiclassical form of the propagators that represent super-operators. 

The following section reviews the semiclassical correspondence between classical 
canonical transformations and quantum unitary operators. The choice of a particular 
representation for the latter is necessarily paired to a specific type of generating function 
for the corresponding classical transformation. In the case of the centre, or the chord 
representations, the respective centre and chord generating functions are subject to 
caustics wherever the canonical transformation can be locally approximated by either 
a reflection, in the case of the centre generating function, or a translation ^Hj- These 
give rise to spurious singularities in the simplest semiclassical approximations that are 
reviewed in section 3, which also presents the corresponding forms of the Hamilton- 
Jacobi equation appropriate to both these generating functions. 

Up to this point, the evolution of the unitary operators results from solving the 
Schroedinger equation. It is only in section 4 that the less studied case of the Heisenberg 
evolution is reviewed with reference to the Weyl propagator and its corresponding 
Hamilton- Jacobi equation. No distinction will be made here with the von Neumann 
equation, appropriate to density operators that are evolved unitarily, because it is 
simply related to the Heisenberg equation by a time inversion, or by changing the 
sign of the Hamiltonian. A direct semiclassical correspondence follows for this type of 
quantum evolution in the case of unitary operators, if we view a fixed classical canonical 
transformation from a moving (canonical) coordinate frame. 

The notion of double phase space is introduced in section 5, within the general 
point of view developed by Amiet and Huguenin |TT| Ej. As well as providing an 
appealing reinterpretation of the previous results, this allows the simplest derivation 
of the Hamilton- Jacobi equation for the chord generating function. A revised study of 
the general evolution of operators follows in section 6. 

Section 7 is dedicated to the semiclassical evolution of translation operators. They 
correspond to planes in double phase space that are analogous to the momentum planes 
in ordinary phase space. Hence, the chord symbol is a 5-function, whereas the Weyl 
(centre) symbol is a plane wave. The latter is the preferred basis in which to view the 
general nonlinear semiclassical evolution of tranaslations into a single WKB-like wave. 
The Weyl symbols of translation operators are propagation kernels for the evolution 
of arbitrary operators that are specified initially by their chord symbol and finally by 
their Weyl symbol. The invariance of the identity operator under Heisenberg evolution 
corresponds to the classical invariance of the zero-chord plane. It is only within this 
plane that the classical double phase space motion generally coincides with the classical 
Liouville flow in ordinary phase space. 

In section 8 we examine the evolution of the conjugate reflection operators. These 
correspond to planes in double phase space which are analogous to the position planes 
in ordinary phase space. Hence, the Weyl symbol is a 5-function, whereas here it is the 
chord symbol that evolves into a single nonlinear WKB-like wave from an initial plane 
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wave. Therefore, it is guaranteed to have no interference and no caustics for a finite 
time. The chord symbol for the reflection operator is the propagation kernel for initial 
Weyl symbols of arbitrary operators to be later specified by their chord symbols. Of 
course, the final Weyl symbol, or Wigner function, which is the subject of parallel work 
[Hj, can be obtained at any time by a Fourier transform. An example is provided of 
evolution driven by a homogeneous cubic Hamiltonian. This is the simplest nontrivial 
case, since the evolution for quadratic Hamiltonians reduces to the linear Liouville flow. 

Our conclusions in the final section emphasise that actually we are dealing with a 
single mixed propagator, in spite of the alternative definitions that correspond to the 
motion of alternative planes in double phase space. Its simple semiclassical form is 
in sharp contrast to the direct centre-centre propagator for Wigner functions, whose 
caustic structure is analyzed in the Appendix. 

2. Correspondence between canonical and unitary transformations 

In this section we recall several known facts about classical-quantum correspondence 
of unitary operators and canonical transformations. This generally only holds within 
a semiclassical approximation, but the reflection and translation operators, which form 
respectively the bases of the Weyl representation and its Fourier transform, belong to 
the special class for which the correspondence is exact. 

Time dependent unitary operators that act on quantum states in Hilbert space 
correspond classically to evolving canonical phase space transformations. In the case 
of motion generated by a constant Hamiltonian operator, H, the continuous group of 
unitary operators, 

U t = exp(-itH/h), (2.1) 

is related to the continuous group of canonical transformations Indeed, if we define 
the points in the 2L-dimensional phase space as x = (pi, ...,Pl, Qi, <?l), we have 
Ct : xq — > Xt, driven by the classical Hamiltonian H(x) according to Hamilton's 
equations, while the Hilbert space vectors evolve linearly: 1^) = U t \ipo). The 
various representations of the unitary operators correspond to different generating 
functions. Semiclassically, the latter determine the phase of the corresponding quantum 
propagators. The Schroedinger equation for U t , e.g. in either of the position, 
momentum, or the Weyl- Wigner representations, corresponds to alternative versions 
of the Hamilton- Jacobi equation (see e.g. [T2|). 

Even an initially single valued generating function may become multiple valued as it 
evolves. Its different branches are then connected along caustics, which lead to spurious 
singularities of the semiclassical approximation for the various representations of U t . 
For this reason, prescriptions for transforming between the different representations is 
an important part of the semiclassical theory ^3]. in the case of the Weyl propagator, 
Ut(x), that is the Weyl- Wigner representation of Ut, the initial generating function 
corresponding to the identity operator /, is just S 4 (x) = 0, so that the appearence of 
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caustics is optimally delayed. Even so, they may eventually appear [TU], occuring where 
C t can be locally approximated as a reflection about a point x, i.e. the transformation 

R x : x_ -y x + = -ar_ + 2x. (2.2) 

It is then appropriate to switch to the complementary phase space representation, 
which, for an arbitrary operator, A, takes the form 

A (0 = ~77T~Tl I dx ex P {'-A A *W). ^ 2 ' 3 ) 



(27Tfi) 

Here we have used the skew product, 

L 

x A x' = (PiQi ~ QiP'i) = J x - x', (2.4) 

n=l 

which also defines the skew symplectic matrix J. In the case of unitary operators, 
the Fourier transform between Z7t(x) and Ut(£) corresponds classically to a Legendre 
transform between St(x) and the new generating function, S t (C,). This describes the 
evolution in terms of uniform translations , 

T% : X- — ► x + = x_ + ^, (2.5) 

by a vector ^, which is a chord in the classical trajectory, hence this is termed the chord 
generating function. 

Corresponding to this, the chord representation [TU], of an operator, A, can 

be defined directly as a superposition of the unitary translation operators 

f 5 = expfVAxV (2.6) 



where x = (p, g) , that is, 



Recalling the definition of the unitary reflection operators as 

and comparing with (|2.3|1 . we see that the Weyl representation, A(x), decomposes an 
arbitrary operator, A, into a superposition of reflection operators [7j: 

A= /" dxA(x) 2 L J R X , (2.9) 

Probably the first to remark on the general structure of translations and reflections 
underlying the Weyl and the chord representations were Grossmann and Huguenin |14j . 

The greater familiarity of the Weyl-Wigner representation is justified by its use for 
the density operator, p, which is represented by the celebrated Wigner function [THj . 
W(x) = p(x)/(27r/i) i , a real (but not necessarily positive) quantum quasiprobability 
in phase space. In the case of the density operator, its Fourier transform can be 
interpreted as a kind of quantum characteristic function, called in [IBJ simply 

the chord function. The evolution of density operators in open systems is constrained 
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Figure 1. The point x-{x) is mapped onto £+(x) by the reflection i? x in the same 
way as it is moved by the given canonical transformation C. It follows that x_(x) is 
a fixed point of the transformation i? x o C, in which case the circuit 74 is defined by 
the composition of both these motions. 



by trace conservation and positivity [3[T7j, but is not limited to the simple form that 
is studied in this paper. It should be noted, however, that the basis operators for 
the Wigner function, namely the reflection operators, i? x , are not themselves positive, 
having degenerate eigenvalues ±1. 

3. The semiclassical propagators 

A given centre generating function, S'c(x), defines the canonical transformation, C : 
x_ — > x+, implicitly through the equation, 

^ = J£c(x), (3.1) 
ax 

with £c(x) = x + — X-, i. e. the chord centred on x [TTJ [TU]- This is subject to the usual 
constraints for generating functions ^H] • It was shown by Marinov ^H] that this centre 
action is just 

5 t (x) =A t (x)-tH(x-), (3.2) 
where A t (x) is the symplectic area 

A t (x)=<fp-dq (3.3) 

for the circuit % that starts along a trajectory segment chosen so that its endpoints 
are centred on x and then is closed by reversing £t(x), as shown in FigQ] Of course, 
for a Hamiltonian system, we have H(x+) = H(x-), but generally H(x) 7^ H(x-). It 
might appear awkward to find the point a?_(x), but this is merely the fixed point of the 
canonical transformation R x o C, which shows this to be a reflection representation of 
classical mechanics ^U], as shown in Figd For short times, the trajectory segment is 
approximately £t(x) itself, so that, like the trajectory, the chord is unique. In this limit, 
the first term in (|3.2j) is negligible, resulting in the explicit approximation 

S t (x) _> o -*if(x). (3.4) 
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Linearizing the canonical transformation in the neighbourhood of the tips of the 
chord centred on x, so as to define the tangent map, 

5x+ = M(x) 5x., (3.5) 

we obtain the Hessian matrix of S(x) as 

^ = -J[l-M][l + M]-\ (3.6) 

which is a Cayley parametrization of M [TT) ITU] . References PUJ E] discuss some of the 
properties and a few examples of the centre generating function. It follows from (J3.6|) 
that the centre caustics for the generating function are the manifolds where M(x) has 
an eigenvalue —1, so that locally C is a reflection for one of the degrees of freedom. 

Let us now consider the evolution equation of a centre generating function, St(x), 
driven by a constant Hamiltonian, H(x). Fixing the initial phase space point x_ from 
which evolves x + (t) = Ct(x_), implies that SSt = —H{xJ)8t, because of the centre 
variational principle [2UJ ED]- Recalling that x_ = x — £/2 together with ()3.1j) . leads 
then to Marinov's version of the Hamilton- Jacobi equation [TUj : 

Several other partial differential equations will be presented in this paper that describe 
the evolution of generating functions, in different contexts and depending on different 
free variables. We here follow the tradition of naming all of these anachronically as 
Hamilton- Jacobi equations. 

These are the basic ingredients for constructing the semiclassical Weyl propagator, 

C/(x) = 2 L | det(l + M)|~3 exp[^" 1 5(x) + W]. (3.8) 

Here, it makes no difference if the unitary operator, U, with Weyl symbol, U(x), 
corresponds to a given canonical transformation C, or a time dependent transformation 
generated by a Hamiltonian H(x). The original derivation j2TJ of U(x) in the latter 
case made no reference to Marinov's generating function, but ()3.8|) follows by inserting 
()3.2|) in Berry's result. An alternative derivation based on the Weyl path integral is 
presented in ^UJ ■ As usual, the semiclassical propagator is exact in the case of quadratic 
Hamiltonians, i. e. linear classical motion. 

For a unitary transformation, Ut, that evolves continuously from the identity, /, 
driven by H, we obtain from ()3.4j) and the fact that Uo(x) = 1 that the phase 9 = for 
small times. This is only altered when a centre caustic of the generating function 5 t (x) 
is reached, so that the denominator of ()3.8j) blows up. But in this event, the eigenvalues 
of M that become —1 must come in pairs, because this is a symplectic matrix (see 
e. g. [ID])- Thus the change of phase can only be a multiple of tt/2, instead of just 
7r/4, which is more usual for caustic traversals. Beyond the passage of the first caustic 
through a given point x, the centre generating function is no longer univalued, so that 
the semiclassical Weyl propagator becomes a superposition of terms like (|3.8j) . one for 
each branch of the generating function. Of course, close to the caustic, the spurious 
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singularity of the semiclassical amplitude points to the need for an improved uniform 
approximation. One way to derive this is to transform back from the chord propagator 
that will now be examined. This is the natural procedure from the double phase space 
point of view that will be discussed in section 5. 

We can now obtain the semiclassical chord propagator by inserting (J3.8)) into (|2.H|) . 
that is 

U(0 = -L, ( dx exp[^(x)-eAx) + ^] (3 g) 

(7rh) L J |det(l + M(x))|5 

and then evaluating this Fourier integral by the method of stationary phase. The 
stationary phase condition simply picks out the centre x(£), such that £(x) given by 
()3.1|) equals the given chord £. The amplitude factor for stationary phase integration is 
then obtained from the determinant of as 

£7(0 = | det(l - M)| - 3 exp[zfi- 1 i S(0 + i</>], (3.10) 

where we identify the chord generating function, £(£), as the Legendre transform of 
S(x). Indeed, it describes the same canonical transformation C by the relation [TU] . 

-|p = -jxc(e), (3 - n) 

which is conjugate to (j3.1j) . As discussed in [TO], we determine £_(£) = x(£) — £/2 
as the fixed point of o C by a construction that is similar to Fig^ Actually, 
there is a conjugate Cayley relation between the matrix M that linearizes the canonical 
transformation and the Hessian matrix of the chord action: 

^ = -J[l + M][l-M]- 1 . (3.12) 

Thus we see that the spurious semiclassical singularity of U(£) in ()3.1U|) takes place at 
the caustic of the chord generating function where det[l — M] = 0. This occurs when the 
canonical transformation in the neighbourhood of the chord tips can be approximated 
by a uniform translation, or the identity. Thus, the chord generating function is singular 
as t — > for a Hamiltonian flow. In contrast, the chord propagator is perfectly regular 
at reflections, which correspond to the caustics of the Weyl propagator. Hence, the 
complementarity of this pair of phase space representations. 

It has been seen that the phase 9 in ()3.8|) is zero for short times. Therefore, the 
phase in ()3.1()j) will be initially —tt/4 times the signature of the Hessian matrix for 
ff(x(£)) see e.g. [I2])- In the case of the simple quadratic Hamiltonians, this signature 
equals 2, for a harmonic oscilator (L = 1) and it equals 0, for the inverted oscillator. The 
explicit form of St(£) for these cases is discussed in QTT . For quadratic Hamiltonians, 
the semiclassical chord propagator is exact just as the Weyl propagator. 



4. Operator evolution 



There is another kind of evolution which may be imposed on a unitary operator, taking 
along with the Weyl and the chord symbols. This is just the usual evolution according 
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to Heisenberg's equation, 
dU'(t) i 



[H',U'(t)}. (4.1) 



dt h 

Here, it is proper to distinguish the external Hamiltonian, H', from the Hamiltonian 
generator, H, which may have been used to define U'(0) = U itself in (12. lj) . Actually, 
(|4.1|) is more commonly employed for the propagation of observables (such as, q, or H) 
or density operators (with a change of sign), rather than unitary operators. Perhaps 
the more intuitive realization is to consider U as an active transformation, whereas H' 
generates by exponentiation, according to (j2.1j) . a continuous group of passive coordinate 
transformations in Hilbert space, V't- Then we obtain 

U'(t) = VL t UV; (4.2) 

as the time dependent unitary transformation that results from the adoption of the 
moving coordinate frame. It is convenient to follow the nomenclature of Osborn and 
Kondratieva [3], so as to distinguish the Heisenberg evolution of U'(t) governed by (|4.1|) 
from the Schroedinger evolution of tl t . Though unusual, the latter is apropriate in 
as much as (j2.1j) is the solution of the Schroedinger equation for the initial condition 
Uq = T. 

This specialization of the Heisenberg evolution to unitary operators is the key 
to a clear classical correspondence. Indeed, we can also view the classical canonical 
transformation, C : X- —>■ x + , in the moving coordinate frame defined by the continuous 
group of canonical transformations, K.' t : x — * x t . The latter is generated by the classical 
Hamiltonian H'{x), that corresponds to H', not to the generator of the transformation 
C. Then the evolving canonical transformation corresponding to U'(t) is just 

C(t) : x'_ -> x' + = K' t o C o K' t (x'_). (4.3) 

Clearly, the evolution of the initial active transformation C depends here on a pair of 
orbits generated by the external Hamiltonian H'(x), the one for x + moving forwards in 
time, while the orbit for x_ moves backwards (or forwards for —H'(xj). This geometry 
is sketched in FigEl The fundamental difference between the classical Heisenberg 
evolution, C'(t), and C t is highlighted by the example where the initial active canonical 
transformation is just the identity: C I. No change of coordinate system can alter this, 
so C(t) = I for all time, which is not generally the case for the canonical transformation, 
Ct, that results from the integration of Hamilton's equations. Of course, the quantum 
operator / is also invariant according to (|4.2jl . in contrast to the nontrivial evolution 
specified by (|2.1|) . 

The semiclassical Heisenberg evolution of operators in the Weyl representation 
was derived independently in 2 J and in [3J. In both cases, the density operator p(t) 
was considered, i. e. the propagation of the Wigner function, but the results can be 
immediately appropriated for the Weyl symbol of U'(t), with a change of sign for the 
time (strictly the von Neumann equation). U'(x.,t) has the same general form as (|3.8|) . 
but with the alternative generating function S"(x, t), corresponding to the canonical 



On the semiclassical evolution of quantum operators 



10 




Figure 2. The classical Heisenberg evolution, C'(t) : x'_ — > x' + , of the initial canonical 
transformation, C : x — » x + , is constructed with the addition of the forward orbit, 
: x + —> x' + , and the backward orbit, K'_ t : x'_ — > x_. 



transformation C'(t) specified by (J4.3|) and the stability matrix M'(x, t) obtained from 
(|3.6|) . The partial differential equation satisfied by S"(x, £) is then a kind of Hamilton- 
Jacobi equation [3], 



95' . „,/ IBS'' 



H' 



X+ 2 J ftT 



0. 



(4.4) 



It will be shown in the following section that the chord generating function for the 
Heisenberg transportation of a canonical transformation satisfies a similar Hamilton- 
Jacobi equation. In the conclusion of the present section this same equation is 
anticipated through the study of the semiclassical limit of the Heisenberg equation in 
the chord representation. So as to represent the commutator in (|4.1|) . we note that (see 

e.g. nm) 



1 



(4.5) 



(2%n) L 

In the case of the Hamiltonian it is better to use the Weyl representation, which is a 
smooth real function. On the other hand, linear phases can be incorporated as a shift 
of the origin [16] : 



A(£)exp(-77A£ 



A{t-r,). 

Hence, the chord commutator can be expressed as 



[H',u\tm) 



(2vrfi) J 



dgdx! 



U , {g,t)e*p[--(Z 



(4.6) 



OAx' .(4. 



Inserting the semiclassical approximation ()3.1()j) for U'(£, t) in this equation, the full 4 L - 
dimensional integral can then be evaluated by stationary phase. First we note that the 
stationary conditions are just that £' = £ and that x' = x(£',t), according to (jH.llj) . i.e. 
the stationary centre that corresponds to the given chord is specified by the classical 
transformation. The Hessian determinant of the phase is unity, so that, within the 
semiclassical approximation, 



H'(3 



2' 



d£ 2' 



U%t). 



(4- 
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In the next section we will verify that the square brackets on the right specify 
the time derivative of S'(£, i). Though it is possible to derive directly the Hamilton- 
Jacobi equation that determines the evolution of the chord generating function for the 
transported canonical transformation ()4.3|) following [2|, the phase space geometry of 
fitting polygons is quite complicated. It is worthwhile instead to invest on the double 
phase space approach followed by [3]. The great advantage is that then chords and 
centres merely define conjugate planes in the enlarged space, so that their relation is 
entirely analogous to the conjugacy of positions and momenta in ordinary phase space. 
In contrast, equation (|4.4jl appears strange indeed when examined in ordinary phase 
space, if it is recalled that ordinary Hamilton- Jacobi equations depend on only half of 
the phase space variables. 

5. Double phase space 

It might seem perverse to double the phase space of classical mechanics, which 
is already a doubling of position space. Nonetheless, we are here concerned with 
representing operators, commonly represented by both bra and ket spaces, so it is not 
surprising that the classical correspondence generally calls for the doubled classical 
space. Observables are deceptively simple in this respect, but, corresponding to unitary 
transformations, there arises an attractive and simple geometrical picture for the 
canonical transformations, C : X- — > x + , defined in the original phase space. Indeed the 
canonical property demands that all closed curves, 7_, be mapped onto closed curves, 
7 + , such that 



Therefore, the definition of the double momentum space, P = (— p_,p + ), and the double 
positions, Q = (g_,g + ), allows us to reinterpret the canonical condition as 



where T = (7_,7 + ). These are arbitrary closed curves on the (2L)-dimensional surface 
defined by the one-to-one function, x + = C(x_), within the (4L)-dimensional double 
phase space X = (P, Q). 

In other words, the action, or symplectic area for any closed curve drawn on the 
surface that defines the canonical transformation, C, in double phase space is zero, i. e. 
canonical transformations are described by Lagrangian surfaces in double phase space. 
This Lagrangian property allows us to define locally a function, 



which is independent of the path followed between Q and Q. In its turn, this generating 
function defines the given Lagrangian surface by the equations 






(5.2) 




(5.3) 



dS c 



Pc(Q), or 



dS c 



dS c 



(5.4) 



dQ 



dq+ 



P+, 



dq— 



= -P-, 
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that determine implicitly the canonical transformation |18j . 

Though the mapping C : a;_ — ► x + = C(a;_) is necessarily univalued, no such 
restriction results on the function Pc{Q), defined by the same (2L)-dimensional surface. 
We cannot define generating functions using x_ or x + as independent variables, because 
these do not lie on Lagrangian planes in the double phase space, i. e. there is 
no constraint that either side of (j5.1|) be necessarily zero. What is allowed and 
often desirable is to apply linear canonical transformations to the double phase space, 
X — > X', which leave invariant the Lagrangian property for any surface, including 
P' = 0. Then we may define a new generating function in the new variables such that 
P' c {Q')=dS' c /dQ'. ^ 

All the commonly used generating functions ^H] are obtained by the application 
of canonical 90° rotations to the single phase spaces, q± — > p±, separately or in 
combination. Obviously, there exist unlimitted other possibilities in double phase space 
[TT] . but we will here be concerned only with the special canonical variables, 

q' = x = ^4^' p ' = y = J ( s + - *-) = J t ( 5 - 5 ) 

Instead of the previous 90° rotations, this transformation to canonized centre and chord 
variables is more like a 45° rotation in double phase space. The plane y = (or £ = 0) 
clearly specifies the identity transformation, I, which corresponds to a Lagrangian plane. 
Actually, all planes y = constant are uniform translations, T^, by the vector — Jy = £, 
whereas each plane defined by a constant x identifies the reflection i? x . Unlike the 
Lagrangian plane (g_,g + ), the planes y = and x = can be considered as phase 
spaces on their own: the space of reflection centres (Weyl space) and the space of 
translation chords. But it must be remembered that these are Lagrangian as far as the 
double phase space action, or symplectic form is concerned. Therefore, the mapping C 
defines implicitly the local function yc( x ) i n terms of the generating function «Sc(x): 
yc(x) = J£ c ( x ) = dSc/dx, which provides a double phase space interpretation for 
the relation between centres and chords in (j3.1j) . Alternatively, the generating function 
Sc(y) can be defined, such that xc(y) = dSc/dy, corresponding to the chord generating 
function ([3. II)) . 

Rather than derive the intricate single phase space geometry for the evolution 
of the chord generating function, it is much simpler to obtain the Heisenberg form 
of the Hamilton- Jacobi equation corresponding to ()4.4|) from the double phase space 
dynamics. First, we note that inserting ()3.1J) into ()4.4j) determines the double phase 
space Hamiltonian as 

H\X) = H'(x, y) = H'(x - Jy/2) - ff'(x + Jy/2) = H'(x+) - H'{xJ). (5.6) 

This should be contrasted to the Schroedinger double phase space Hamiltonian obtained 
from (|3.7|) . which is simply 

M(X) = H{x, y) = if (x - Jy/2) = H(x+). (5.7) 

In each case IH(X) generates the Hamiltonian flow in double phase so as to evolve 
any Lagrangian surface acording to the respective Hamilton- Jacobi equation. The 



On the semiclassical evolution of quantum operators 



13 



Heisenberg double Hamiltonian (J5.6j) determines the phase in Marinov's path integral 
for the Wigner function [22J. 

It must be stressed that the motion generated by the double Hamiltonian, M'(X) 
in (|5.6|h is purely classical, albeit in double phase space. This is in general quite different 
from the Liouville flow in single phase space. Generally, it is only within the invariant 
plane, y = 0, or £ = 0, that the flow is Liouvillian, as will be shown in section 7. 

Generally, the double phase space motion generated by IH'i^x, y) in ()5.6|) depends 
on both the initial centre, x and on y . Indeed, double phase space can be considered 
as an incorporation of the dynamics of pairs of orbits, starting from x_o and x + o, into a 
single Hamiltonian scheme. Thus, it is only within the 'horizontal' zero chord plane of 
double phase space that the single and double phase space motions always coincide, as 
will be further discussed in the next section. However, in the special case of quadratic 
Hamiltonians, H'(x) = —a A x + x ■ Bx, where B is an orthogonal matrix, the motion 
of chords and centres becomes independent. This follows simply from (|5.6|) . so that 

H\x,y) = ay-2x-BJy, (5.8) 

and Hamilton's equations for the conjugate variables x and y become 

x = — — = -a - 2JBx (5.9) 
dy 

and 

dm' 

y = — = -2BJy. (5.10) 
ax 

In the case of the centres, this is exactly the same as x in single phase space and merely 
reflects the fact that the centre between two points satisfies the same linear equation 
as each of them individually. This fact is the basis of the classical propagation of the 
Wigner function and Weyl symbols for quadratic Hamiltonians. What is not so familiar 
is that the chord function and chord symbols also propagate classically [23] ■ In the 
case of a homogeneous quadratic form (a = 0), Hamilton's equations for £ = Jy and 
for x again coincide according to [23J, but the chords are insensitive to the linear part 
a A x. The latter generates a translation, which appears in the chord symbol as the 
phase factor in (|4.fi|) . 

Coordenatizing locally the moving surface driven by 1H'(X) as y = y'(x, t), the 
solution is the generating function 

S"(x,*) = / X y / (x,t)-dx. (5.11) 

The conjugate action is defined via the Legendre transform, 

S'(y,t)=x!-y-S , (x!,t), (5.12) 

in the usual way, i. e. x'(y,t) is specified by the condition that ()5.12|) is stationary 
with respect to x'. Taking now the full derivative of ()5.12)1 with respect to the time and 
using the fact that dS'(y,t)/dy = x'(y,t), then leads to 

^) = _^) = « V(y , t ), y) . (5 ,3) 
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Reintroducing (|3.11|) into (|5.13|) and again recalling that £ = — Jy, leads to the chord 
evolution equation corresponding to the Heisenberg evolution as 

Of course, we can also obtain the partial differential Hamilton- Jacobi equation for the 
Schroedinger action S t (x) by exactly the same procedure, so that 

6. Operator evolution revisited 

To apply the above theory to the evolution of quantum operators it should be recalled 
that they form a Hilbert space of Hilbert-Schmidt operators with the scalar product 

((A\B)) = tr$B, (6.1) 

defined in terms of the adjoint operator, A'. When the trace is defined, the Heisenberg 
(or von Neumann) evolution can be considered as the action of a unitary super- operator, 
since it preserves the scalar product. General super-operators always preserve the trace 
of the self-adjoint density operator, trp, but not necessarily trp 2 

Each foliation of double phase space by parallel Lagrangian planes corresponds to 
a possible operator representation. Perhaps the most common representation relies on 
the position projection operators, ((Q\ = |g_)(g + |, so that 

((Q\A)) = (q + \A\q.)=tr\q.){q + \A, (6.2) 

where the Lagrangian planes are just Q = (g_,g + ) = constant. From this one can 
switch to momentum, or various mixed representations through Fourier transformations, 
corresponding to 90° rotations in double phase space. The Weyl representation, based 
on the self-adjoint operator, i? x , then corresponds to the double phase space rotation , 
so that for Q' = x, 

((Q'\A)) = 2 L trR x A = A(x), (6.3) 

and the Lagrangian basis in double phase space are the reflection planes, x = constant. 
The Fourier transformation ()2.8|) then brings in the translation operator, whose adjoint 
is T_£. This is represented in double phase space by the new Lagrangian plane, 
P' = y = J£ = constant, so that 

((P'\A)} =tr f^A = A(0. (6.4) 

In each case the representation in terms of a set of Lagrangian planes, Q' , is 
complementary to the conjugate representation in terms of P' , which is obtained by 
a Fourier transform. Thus, Heisenberg's uncertainty principle manifests itself in double 
phase space. Further discussion is presented in reference 
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The suggestive use of a Dirac notation for operators, in direct analogy to that 
commonly reserved for states in Hilbert space, is supported by the manner in which the 
evolution of unitary operators is semiclassically related to the evolution of Lagrangian 
surfaces, in double phase space. In the case of the chord propagator, combining (|5.14|) 
with ()4.8j) we verify that 

[H\U\t)m = ~^(U), (6.5) 

to lowest order in Planck's constant. In the following sections we study the evolution of 
the particular unitary operators, translations and reflections, which are the basis of the 
present theory. 

The preceding theory is concerned with the evolution of unitary operators that 
always correspond to Lagrangian surfaces in double phase space. The next section 
focuses on evolved unitary translation operators which allow us to propagate arbitrary 
operators through their chord symbols. But first, it should be noted that the present 
semiclassical theory can also be applied directly to other operators that also correspond 
classically to Lagrangian surfaces in double phase space. 

The most important case is that of projectors, or pure state density operators, 
|V>)(0|- Indeed, it is quite usual for the semiclassical description of a state \ip) to be 
based on a Lagrangian surface in single phase space, a torus for a bounded state, e.g. a 
closed curve in the case that L — 1. Then the projector onto this state is supported by 
the product manifold with doubled dimension, e.g. a two-dimensional torus if L — 1. 
Another interesting case is that of a dyadic operator, \ip){ip'\, which is useful to describe 
transitions. In the case that and l^') are eigenstates of the same observable, the Weyl 
representation of this operator is known as a Moyal fuction, or cross- Wigner function 
[2111211 • Its semiclassical form in single phase space has been previously studied [28J, 
relying on chords (and their centres) that connect the pair of distinct single tori, but the 
double phase space picture is simpler because a single Lagrangian surface is involved as 
in the case of the projectors. 

In all these cases, caustics cannot generally be avoided. Taking the product of 
the torus in x_ with the torus in x + , results in a surface with the same dimension 
as the single phase space in which each torus is embedded, but which projects down 
singularly as the single torus into either of x± . But, as we have seen, x± are not 
Lagrangian surfaces in double phase space with which to view caustics of the double 
torus in the sense that they cannot supply the variables for a generating function. The 
£ = plane within double phase space, in which the Wigner function is defined, is 
Lagrangian and the caustic that arises along the torus in the single space picture 
can now be reinterpreted as the fold caustic for the projection of the double torus onto 
this particular coordinate plane. The section of the double torus with the plane £ = 
is identical to the single space torus, just as the sections with x_ and x + . If L > 1, 
the Wigner caustic has a higher dimension than the single torus, though the latter is 
included as a higher singularity j30j- In all cases, the Wigner caustic results from the 
projection of the double torus onto the x plane. Further application of semiclassical 
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methods in double phase space are presented in 

The projection of the double torus onto the chord plane, x = 0, is even more 
singular. Indeed, the origin of this plane, £ = 0, is the image of the entire single torus 
The approximation of the chord function in this region is studied in ^H], as well as the 
caustic for maximal chords of the single torus. In any case, it is clear that a simple 
Fourier switch between the Weyl and the chord representations does not deliver us from 
caustics and the problem of their evolution. The alternative that is pursued in sections 
7 and 8 is to propagate either the translation, or the reflection basis operators. 

7. Propagating translations 

It was reviewed in section 2 that the chord representation may be considered as the 
decomposition of an arbitrary operator into a superposition of unitary translations. 
Hence, the general evolution of operators in this representation can be reduced to the 
propagation of this special class of operators, each of which corresponds to a basis 
plane y = constant in double phase space. These can be pictured as horizontal, if we 
interpret y as the double momentum, P', and the y = £ = plane corresponds to the 
identity transformation. Each of its points results from the intersection with a vertical, 
Q' — x = constant, plane that corresponds to a reflection, being that the plane x = 
corresponds to the operator for reflection through the origin of ordinary single phase 
space, which is also known as the parity operator. 

It is immediately obvious from the form of the double Hamiltonian for Heisenberg 
propagation (|5.fijl that the y = plane is always invariant. This merely expresses 
the fact that the corresponding canonical transformation remains the identity for all 
time, just as the quantum identity operator, /, is invariant for arbitrary Heisenberg 
propagation. Therefore, classical-quantum correspondence becomes exact in the limit 
of small chords. This confirms the interpretation in JH] of the long chords as responsible 
for the essentially quantum long range correlations of the density operator. Within the 
purely quantum domain, the combination of the unitarity of the Heisenberg super- 
operators, which preserves the trace of products of operators (jfi.lj) . with the invariance 
of J, results in the invariance of the trace itself. This property must be demanded of 
the nonunitary super-operators for open systems, which are no longer generated by the 
Heisenberg equation jHj. 

The flow within the invariant classical plane results from the expansion of the double 
phase space Hamiltonian, M' (x, y) in (|5.fij) . for a fixed centre x, 



Hence, Hamilton's equations in double phase space for the centres, x, within the classical 
invariant plane are identical to the equations of motion generated by the single phase 
space Hamiltonian, H'(x). 

General translations, C = T^, are not invariant for arbitrary Heisenberg 
propagation. The initial Lagrangian surface that will be evolved by C'(t) = T'At) is 




dH' 
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a horizontal plane, so that analogy with single phase space indicates that the chord 
symbol is 

r € (0 = *(e-0> ( 7 - 2 ) 

whereas the Weyl representation is the plane wave [TUj . 



7Hx) = exp(--xA£). (7.3) 

The chord symbol is not in a simple semiclassical form, because the horizontal planes 
project singularly onto the vertical axis, but (|7.3|) is in the form ()3.8|) with the classical 
action 5^(x) = — x A £ and the tangent matrix M^, defined by (|3.5j) is twice the unit 
matrix. Inserting this linear form of the generating function in (|3.1|) results in the same 
chord £ being placed on all the centres x. This need no longer hold after a classical 
Heisenberg evolution, so that the evolved actions ^(x', t) obtained from the Hamilton- 
Jacobi equation (|4.4|) will generally have second and higher order terms and M^(x',t) 
will no longer be proportional to the unit matrix. Unless the evolution proceeds to the 
possible production of a vertical fold, i.e. a centre caustic, the phase 6 in ()3.8|) will be 
the same as in (|7.3|) . Thus the general semiclassical form of the distorted translation in 
the centre representation is 

T 5 '(x',t) =2 L |det(l + M^(x , ,t))|-^exp[^- 1 ^(x',t)]. (7.4) 

The identification of the Hermitian conjugate of T^(t) with T"_^(t), i.e. the inverse 
operator, establishes that S'_g(x,t) = —Sg(x,t). 

An important constraint on the classical Heisenberg evolution of the plane that 
describes a finite translation is that it cannot touch, or intersect the invariant identity 
plane, y = 0. This means that the distorted translation never develops a fixed point. 
Indeed, the number of fixed points is invariant throughout the evolution. Furthermore, 
the linearization of the monodromy matrix near each fixed point is also invariant. Thus, 
summing over the fixed points we obtain a semiclassical invariance of Tabor's version 
[3*Tj of Gutzwiller's trace formula [32] , once the quantum invariance of the trace of any 
operator is recalled. 

If we now adopt the Weyl representation for the Heisenberg evolution A'(t) of an 
arbitrary operator A, ()2.7|) leads to 

il 

(2vrfi) J 

Hence, the distorted translation, T^(x!,t), is identified as the mixed chord-centre 
propagator. In view of the previous discussion, the simple semiclassical approximation 
()7.4j) then supplies an initially caustic free approximation for the Weyl symbol (or 
Wigner function) that evolves from a given chord symbol. This is not the case of 
the chord-chord propagator that follows from the evolution of the chord symbol of A'(t) 
given by ([2.7)1 : 

A'tf,t) = f 1 0y;A{t)rtf,t). (7.6) 



/ct^)*.*)- (7-5) 
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Even though the invariance of trA'(t), implies that A'(0,t) = A' (0,0) for all t, 
this propagator evolves from the 5-function (|7.2|) into a semiclassical form that must 
negotiate caustics through Airy functions, or other higher diffraction catastrophes 
[33J[T2|. Using the product rules for chord symbols [TU] it is possible to express the chord- 
chord propagator directly in terms of the chord symbols for the elementary evolution 
operators V{ in the full Heisenberg propagation (|4.2|) : 

rtf, t)=J d V V t ' L + VI (rj - * exp (U) A *±?j .(7.7) 

It is interesting to note that this exact expression has the same form as the 
transformation that defines the chord function x(Q f° r a pure state density operator 
from the wave function (q\ip) [TO] . Stationary phase evaluation of (|7.7|) far from caustics 
leads to a superposition of terms of the general semiclassical form (|3.1U|) . 

As noted at the end of section 6, we are dealing with a general semiclassical theory 
for the evolution of operators, each of which corresponds to a Lagrangian surface in 
double phase space. The translation operators are an important particular case, but if 
the operator A that is propagated is related to a particular Lagrangian surface of its 
own, one may as well evolve it directly in the Weyl representation using (|3.8|) . instead of 
integrating the evolved basis operators. Indeed, far from a caustic, the stationary phase 
evaluation of (|7.5|) . or ()7.6j) will give the same result as (|3.8jl . However, specially for 
pure state density operators, the whole region of small chords, that is identified with the 
near-classical part of the density operator [TH], lies in the neighbourhood of a caustic 
of the chord function as well as a caustic of the Wigner function. Furthermore, the 
phenomenon of decoherence, i.e. the evolution of the density operator in an open system 
(into a mixed state) [T7], leads to a cut-off for the contribution of the longer chords [23J. 
It is thus important to determine the semiclassical evolution of the neighbourhood of 
the invariant plane y = (or £ = 0) in double phase space. The difficulty is that 
it is necessary to keep track of the (small) chords throughout the evolution, whereas 
the chord-centre propagator in this section depends on the initial chord, but is later 
only described by the centres. The situation is clarified in the following section by the 
alternative definition of the mixed centre-chord propagator. 

8. Propagating reflections 

The Heisenberg evolution of arbitrary operators, A, according to (|4.2|) . is given by (|2.9|) 
as a superposition of evolved reflection operators. (Though density operators propagate 
backwards in time.) Taking the Weyl symbol of the evolved operator A'(i), then leads 
to 

A'(x',t) = J (ixA(x) 2 L R' x (x',t), (8.1) 

where we introduce in the integrand the Weyl symbol for R' x (t), the quantum reflection 
through x seen in a moving coordinate frame. Thus, (|8.1|) identifies 2 L R' x _(x f ,t) as 
the propagator for Weyl symbols and Wigner functions (with backward propagation in 
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the latter case). This propagator can be expressed exactly as a path integral [22] and 
can also be obtained exactly [2 from the Weyl symbols of the operators V( that are 
responsible for the Heisenberg evolution ()4.2|) . using the Weyl product rules [27)] : 

2X(x', t) = Jdx VI (^±^ + x) VI - x) * exp (|x A (x' - x)) . (8.2) 

This formula is the companion of ()7.7)1 and it has the same structure as the transform 
that defines the Wigner function for a pure state in terms of the corresponding wave 
function [To] . 

Initially, for t = 0, (J8.1)) reduces to the known result that 

ik(x') = 2- L 5(x-x'), (8.3) 

i.e. the Weyl symbol for the reflection through the point x is just a 5-function [7j. 
This is analogous to the fact that position states in single phase space are represented 
by (^-functions in the position representation, because we associate reflection centres to 
positions, Q' = x, and canonized translation chords, P' = y = J£, to momenta. Thus, 
it is no surprise that in the conjugate chord representation we have 

J R x (0=expQxAe), (8.4) 

i.e. the reflection operator is represented by a plane wave JU]- Note that the 
double phase space relation between reflection operators, viewed in the translation 
(chord) representation ()8.4jl . and translation operators in the Weyl (reflection centre) 
representation is in perfect analogy to that between the relation for position eigenstates 
viewed in the momentum representation and momentum eigenstates in the position 
representation. 

In the special case where the Hamiltonian, H', is quadratic, so that the 
corresponding classical transformations are linear, the propagator for Weyl symbols 
remains a Liouvillian 5-function for all time, 

R' x (x,t) = R x( t)(x) = 2- L 5(x(t) - x'), (8.5) 

where x(£) is the classical trajectory for x. The double phase space point of view 
shows how special this situation is: it is a canonical transformation that takes vertical 
planes into vertical planes. The Appendix discusses how any nonlinearity immediately 
introduces caustics into this centre representation, because the transformed surface 
representing the canonical transformation, R' x (t), must have a vertical tangent at y = 0, 
as shown in Fig. 3(b). The nonlinear distortion of a reflection in the original single phase 
space is sketched in Fig. 3 (a). 

In contrast, using the chord representation for these propagated reflections, R'(t), 
the Heisenberg evolution merely distorts the plane wave (|8.4|) . which can only develop 
a caustic after a finite time. Even so, the invariance of the number of fixed points for 
classical Heisenberg evolution commented in the previous section prevents the evolved 
vertical plane in double phase space from intersecting, or even touching the invariant 
plane y = a second time. Comparing ()8.4j) with (j3.10|) shows that the exact R' K {£,', 0) 
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Figure 3. (a) A nonlinear classical Heisenberg evolution evolution maps the original 
reflection centre, x — > x'(i), but the linear reflection about this point is only a local 
linear approximation to the evolved transformation, R' x (t): The centres, x'(£'), of finite 
chords do not generally coincide with x'(t). However, for all chords, x'(£') = x'(— £'). 
(b) In double phase space, a linear reflection is represented by a vertical plane. This is 
curved by a nonlinear evolution, but the tangent is still vertical at the identity plane, 

y = o. 



is already in its semiclassical form, with S' x (£',0) = xAf. The classical Heisenberg 
evolution ([4.3)1 . i.e. the Hamilton- Jacobi equation (|5.14|) . takes this into 

S;(£',*)=x(i)A£' + 0(£' 3 ), (8.6) 

which is still an approximate linear reflection about x(t) for small chords. All even 
powers in the components of £ are missing in S' x (£',t), because the involution property 
of the reflections, that is R\ = I, the identity, remains invariant for coordinate changes; 
hence x(-f,t) = x(f,t)). 

There is a simple semiclassical form for the mixed centre- chord propagator that 
defines the chord symbol A'(£', t) for the operator A'(t), which has evolved according to 
the Heisenberg equation from the operator A with Weyl symbol A(x): 

A'(?,t) = /dxA(x)2X(£',t). (8.7) 

The mixed propagator has the semiclassical form defined by a single evolving classical 
generating function, S' x (£',t), which is a special case of (j3.10|) : 

K(e,t) = |det(l-M;(£',t))|-iexp[^ 1 ^(e',t)], (8.8) 

where the monodromy matrix M' for the transformation from x — £/2 to x + £/2 
is obtained from the evolving action in (|3.12j) . The familly resemblance between the 
semiclassical expressions f|8.8|) for R' x (£',t) and (|7.4|) for T^(x',t) can be pushed all the 
way by the identity, 

7j(x,t) = trf[{t)2 L R x , = tr2 L R' x (-t)f^ = R^,(£, -t), (8.9) 

where (Q is used for f^t) and R' x {-t). It follows that S' x (£',t) = S'^fa-t) and 
hence both semiclassical propagators will be free of caustics for the same interval of 
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time. The identity (|8.9|) expresses the fact that Heisenberg evolution defines unitary 
super-operators, so that a complex conjugate kernel is obtained by reversing the time. 
This unitarity is also apparent in f)8.2j) . 

In the limit of small chords the amplitude in (j8.8|) tends to one, just as in (18. 4J) . In 
fact, it is now legitimate to investigate this limit and thus keep only the lowest term in 
()8.6|) . since £' is the free variable. In this region even the distorted reflexion operator 
takes the simple form 



Here, x(t) is the Liouville orbit of x(0) = x in single phase space, because we are 
restricting the double phase space motion to the invariant plane y = 0, where the double 
Hamiltonian is just ()7.1|) . Thus, the mixed centre-chord propagation (|8.7j) becomes 
approximately 



The small chord approximation is thus a moving Fourier transform of the initial Weyl 
symbol. At first sight, this may appear to be equivalent to generalizing the ^-function 
centre-centre propagator ()8.5|) for nonquadratic Hamiltonians, but this is not so. Such 
a stronger approximation is equivalent to extrapolating (|8.11j) for all chords, which is 
not necessary in the chord representation. 

Consider, for instance, the evolution of the pure state density operators studied in 
[TH] . Even though such operators may be associated to a double phase space Lagrangian 
surface, given locally by the functions £ = £(x), or x = x(£), both the chord and the 
centre representations have caustics in the part of this surface for which £ is small. 
Reference ^H] developed a small chord approximation which could be extended so as to 
overlap with the semiclassical approximation valid for long chords. Now we find that 
these different regions can be propagated separately. For the long chords that carry 
information concerning long range quantum correlations, the simplest is to procede 
with a direct semiclassical evolution, based on the classical propagation of the generating 
function, S" (£,£), for the Lagrangian surface corresponding to the density operator. As 
for the small chords in the near classical region, the best alternative is to use a chord 
basis, so that the evolution is pictured as a moving Fourier transform from the initial 
Wigner function (jS.llj) . 

8.1. A simple cubic case 

The case of a quadratic Hamiltonian is exactly solvable, as mentioned previously. An 
exact result in the nonquadratic case is also available for the Hamiltonian 



We can indeed make the full development of the double phase space Heisenberg 
Hamiltonian (|5.fij) . 




(8.10) 




(8.11) 



H\x) = ap 3 . 



(8.12) 



H'{X) = H'{x + i/2) - if (x - i/2) = a (?>p% + ^ 



(8.13) 
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However, it must be remembered that the relevant classical motion occurs in the 
double phase space X = (x, y = J£), so that the conjugate pair of variables are 
(x p = p, y p = — £ g ) and (x 9 = q, y q — £ p ). Thus, Hamilton's equations for M\X) result 



in constant values for p and £ p , whereas 
= g + 3a(p 2 



~£ 2 )t 



(8.14) 



and £ q (t) = £g + Qap^ p t. Note that the motion within the invariant plane, £ = 0, is just 
the classical motion in single phase space: 

q(t) = q + 3apH. (8.15) 

In this simple case, we can verify that the chord symbol for the reflection through 
the point x, 

at 



satisfies 



dt 
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(8.16) 



(8.17) 



by performing the integrations in the exact expression for the commutator (|4.7J1 . 
Furthermore, the phase 



at 



s' x (£,t)= P r- q (t)e p + -C 



4 -p 



l.tt 



in ()8.16|) satisfies the Hamilton- Jacobi equation (|5.14|) . so that the semiclassical form of 
the mixed propagator (J8.8j) is exact in this example, given that the amplitude is constant. 
Nonetheless, the motion x' (£',£) obtained from the derivative of S' x (£',t) according to 
(ETTTJ) . i.e. (jHH), is not Liouvillian unless £' = 0. 

Inserting (|8.1fij) into the general evolution formula (|8.7jl results in 

A'(£,t) = |dxA(x)2exp [/ / " ' " f 



(8.19) 



Thus, the Fourier kernel that initially transforms from the Weyl representation to the 
chord symbol becomes non-Liouvillian, as well as nonlinear, and the cubic term in the 
phase cannot generaly be neglected. 

Let us now consider the propagation of the Wigner function for the position 
projector |go)(<?o|> which is just Wo(x) = S(q — q ). This corresponds classically to 
the evolution of the vertical straight line which is bent into a parabola under the action 
of the Hamiltonian (|8.12|) . According to ()8.19|) the evolving chord function will be 



Xo(C> *) = / d Pdg 5(g - q Q ) 2 exp 
which can be immediately integrated to yield 
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This is the exact chord function for the parabola, though it is in its semiclassical form 
[To] , so the exponent is just the chord action. Unlike the corresponding Wigner function, 
there is no interference, because a parabola translated by £' intersects the original 
parabola at a single point. In contrast, a parabola reflected through x' intersects twice, 
if x' is in the concave region, or not at all. The corresponding Wigner function, Wq(x', t), 
is an Airy function, as can be verified by directly integrating the Fourier transform that 
relates these representations, using (|8.21j) . The parabola itself is the Wigner caustic, 
whereas the chord function has a nongeneric caustic for £' p = 0, which corresponds to 
the translated parabola approaching the original parabola at infinity, for any horizontal 
translation. 

The small chord approximation (j8.11|) in this case merely misses the (at/4)f£ 3 part 
in the exact expression (j8.21j) for Xo(£'>^)- However, it is a good approximation for 
small chords even along the chord caustic. It is only if this locally valid expression 
is extrapolated for all chords that its Fourier transform yields the crude Liouvillian 
approximation to the Airy function evolution: 



9. Conclusion 

Previous semiclassical approximations for the evolution of the Wigner function, or 
other oscillatory Weyl symbols have not dealt with the near-classical region where the 
conjugate chords are small. However, this will be the dominant region in the integral 



We have shown that the simplest way to deal with the propagation in the near-classical 
region is to resort to the mixed propagators between the chord symbol, A(£), and 
the centre, or Weyl representation. These propagation kernels are defined by the 
Heisenberg evolution of different operators, 2 L R'^(( i ' ,t) and Tg(x',t), namely reflections 
and translations, but they are identified by the relation (|8.9)1 . 

Though it is possible to dispense with the construction of a double space, the 
conceptual clarification also simplifies the calculation. From this point of view, the 
transition between centres and chords lies in strict analogy to that which relates the 
conjugate positions and momenta in ordinary phase space. From this point of view, we 
have merely applied the general Maslov method for negociating caustics to double 
phase space. It turns out that caustics in the evolution kernel are not avoided by 
viewing the entire evolution in neither the Weyl, nor the chord representation on their 
own, but we are guaranteed a simple semiclassical form by alternating between them. Of 
course, the more intricate chord-chord, or centre-centre propagations lie a mere Fourier 
transform away. 

Whether considered in double, or in single phase space, the propagation of unitary 
operators, such as the translations, or reflections which are the bases of the chord 





that determines the expected value of any observable, A: 





On the semiclassical evolution of quantum operators 



24 



and the Weyl representations, have a clear classical analog. This classical Heisenberg 
evolution continuously distorts a fixed canonical transformation through an evolving 
change of coordinates in phase space. Viewed in double phase space, this is an evolving 
Lagrangian surface. Once this quantum-classical correspondence is established, we are 
free to also propagate other Lagrangian surfaces in the doubled space, which, for instance 
correspond to projectors, or dyadic transition operators. In all cases, the motion is driven 
in double phase space by a classical Heisenberg Hamiltonian, simply related to the single 
Hamiltonian by (|5.fi|) . Each trajectory in double phase space corresponds to a pair of 
trajectories in single space, which emmanate from the pair of points x± = x±£/2. Only 
in the limiting case where £ —>■ does the motion of the centre x(£) coincide with the 
single phase space trajectories, x(t), of the Liouville flow. 

Translation chords and reflection centres are conjugate classical variables that define 
transverse foliations of double phase space. They form the bases for complementary 
representations of quantum operators. The extension of Heisenberg's uncertainty 
principle to double phase space, discussed in sections 5 and 6, prevents us from defining 
centres and chords simultaneously within the quantum theory. This fact lies at the 
root of the dificulty of discussing the near-classical propagation of Wigner functions, 
because this region is defined as the region of small chords, whereas the Wigner function 
is constructed in the conjugate centre basis. In references [3 E] this problem was partly 
circumvented by analyzing directly the evolution of given operators, semiclassically 
linked to particular Lagrangian surfaces, which locally tie a single chord to each centre. 
However, the classical region of small chords for all such surfaces are caustics in both the 
chord and the Weyl representation, so that improved uniform approximations become 
necessary. Here, we have adopted the alternative of developing mixed propagators which 
bypass the uncertainty principle because the chord and the centre are each specified 
at a different instant. These propagators are privledged representations of the super- 
operators that act on the space of linear operators of Hilbert space. 
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Appendix A. Centre-centre propagation 

The propagator, 2 L i?^(x',t), that takes a Weyl symbol of an operator into a new 
Weyl symbol and thus evolves arbitrary Wigner functions corresponds classically to 
a Lagrangian surface in double phase space, y^(x', i), that has evolved from the initial 
vertical plane, x = constant. This is strictly analogous to the position propagator, 
(q + \Ut\q-), in single phase space, which corresponds to the classical evolution of the 
vertical plane, q = g_. In both cases the initial vertical plane is an extreme nongeneric 
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caustic, but the semiclassical position propagator is usually well behaved for a finite 
t > 0. It might not then be immediately obvious why caustics are inevitable in the 
double space evolution. In other words, why is the evolving Lagrangian surface allways 
tangent to a vertical plane at y = 0, as shown in Fig.3(b)? 

The essential difference lies in the driving Hamiltonians. The typical form of 
the Hamiltonian in single phase space corresponding to the Schroedinger equation is 
p 2 /2 + V(q), so that q = p. This tilts the initial plane, q = g_, and so breaks the 
verticallity. Indeed, to first order in time, the Lagrangian surface is p = [q — qJ) /t. 
In contrast, the double phase space motion is driven by the Heisenberg Hamiltonian, 
y), defined by (|5.fij) for an arbitrary single space Hamiltonian H'(x). Its expansion 
to lowest order in y is given by (|7.1jl . which leads to an evolution of the x coordinate 
that is independent of y. The third order terms in y contribute quadratic terms to the 
classical equations of motion which break the strict verticality, but the vertical tangent 
remains at y = 0. 

In the example at the end of section 8, the Fourier transform of (18.16)1 leads 
immediately to an Airy function for the centre-centre propagator and generically we 
can expect uniform approximations based on Airy functions in the case of a single 
degree of freedom. This was previously derived for special quantum maps |5%| . For 
higher dimensions and for nongeneric cases, higher diffraction catastrophes come into 
play |3Sj- In contrast, the mixed propagator remains caustic free in all these cases. 
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